Observing complete gravitational wave signals from dynamical capture binaries 
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We assess the detectability of the gravitational wave signals from highly eccentric compact bina- 
ries. We use a simple model for the inspiral, merger, and ringdown of these systems. The model 
is based on mapping the binary to an effective single black hole system described by a Kerr met- 
ric, thereby including certain relativistic effects such as zoom-whirl-type behavior. The resultant 
geodesies source quadrupolar radiation and, in turn, are evolved under its dissipative effects. At the 
light ring, we attach a merger model that was previously developed for quasicircular mergers but 
also performs well for eccentric mergers with little modification. We apply this model to determine 
the detectability of these sources for initial, Enhanced, and Advanced LIGO across the parameter 
space of nonspinning close capture compact binaries. We conclude that, should these systems exist 
in nature, the vast majority will be missed by conventional burst searches or by quasicircular wave- 
form templates in the advanced detector era. Other methods, such as eccentric templates or, more 
practically, a stacked excess power search, must be developed to avoid losing these sources. These 
systems would also have been missed frequently in the initial LIGO data analysis. Thus, previous 
null coincidence results with detected gamma-ray bursts cannot exclude the possibility of coincident 
gravitational wave signals from eccentric binaries. 



PACS numbers: 95.30.Sf, 97.60.Lf 
I. INTRODUCTION 



In dense stellar regions, such as galactic nuclei or glob- 
ular clusters, individual black holes (BHs) or neutrons 
stars (NSs) can become gravitationally bound as en- 
ergy is lost to gravitational radiation during a close pas- 
sage. These dynamically captured pairs may be addi- 
tional sources for gravitational wave (GW) detectors, as 
well as sources of electromagnetic (EM) transients such as 
short gamma-ray bursts (SGRBs). Eccentric pairs will be 
distinguishable from quasicircular coeval binaries, which 
are born in a bound system and have had time to circu- 
larize before reaching the sensitive bandwidths of ground- 
based GW observatories such as LIGO [J, VIRGO 0, 
GEO600 0, and KAGRA(LCGT) 0. 

The primary purpose of this paper is to study the de- 
tectability of sources that retain eccentricity while in the 
LIGO band (for simplicity, we only employ LIGO sen- 
sitivity curves). Before getting into the details of our 
model and results, we briefly review current event rate 
estimates, high-eccentricity population fractions, possi- 
ble EM counterparts, and the GW detectability of dy- 
namically captured compact objects. 
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A. Event rates 

Galactic nuclei are a promising setting for the forma- 
tion of dynamical capture binaries. Mass segregation 
around a central massive BH can lead to large densi- 
ties of stellar mass BHs and stars. For example, the 
Fokker-Planck model used in [5] suggests that our galac- 
tic nucleus should have ~ 2000 BHs and ~ 400 NSs in 
the central 0.1 pc. In 0(7], the event rate for the for- 
mation of BH-BH binaries from GW capture in this set- 
ting was estimated to be roughly between 0.01 and 1.0 
yr _1 Gpc~ 3 , with corresponding Advanced LIGO detec- 
tion rates of w 1 — 10 2 yr^ 1 . This rate assumes that the 
number density n of BHs in galactic nuclei has a scatter 
with (n 2 )/(n) 2 = 30. Assuming no scatter would reduce 
the above rate by a factor of 30. This also assumes a num- 
ber density of contributing galaxies of 0.05 Mpc~ 3 ; i.e. 
it includes all galaxies as contributing roughly equally. 
Lower mass galaxies are not as well understood, though 
if a significant number of them have total cluster mass 
fractions above the 2.5% used in the aforementioned cal- 
culation, this rate would increase. Other unaccounted 
for effects, such as steeper profiles from light-dominated 
mass functions [8], could also potentially increase this 
rate. The formation of BH-NS binaries is estimated to 
be ~ 1% of this rate [5J. 

Dynamical capture binaries may also form in glob- 
ular clusters (GCs) that undergo core collapse [9j fTU] . 
In pj] binary formation through tidal capture was stud- 
ied. Using M15 as a prototypical GC, it was calculated 
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that the NS-NS tidal capture rate would peak at ~ 50 
yr" 1 Gpc~ 3 at z = 0.7 (falling to ~ 30 yr _1 Gpc~ 3 by 
z = 0) for their default model of core collapse. They 
also provide a scaling to BH-NS and BH-BH mergers 
which (assuming Mbh = 4.5M© and a relative fraction 
of BHs to NSs /bh//ns ~ 0.28) gives rates that peak at 
- 70 yr^Gpc" 3 and - 20 yr^Gpc- 3 for BH-NS and 
BH-BH mergers, respectively. This scaling does not in- 
clude complications due to BH ejection [T2TH5] . Also, 
these calculations do not include the likely reduction in 
compact object (CO) populations within the GC due to 
natal kicks. In [T7] it was found that including a 5% 
NS retention fraction when fitting simulation results to 
observations of M15, and assuming no central BH, re- 
duced the estimated number of NSs in the inner 0.2 pc 
by ~ 1/4 compared to a similar study that did not in- 
clude natal kicks [18]. The calculated NS-NS merger rate 
is quite sensitive to the fraction / of NSs in the core, 
scaling as ~ f 2 , which means the aforementioned rates 
could be too large by an order of magnitude if retention 
rates are this low. However, observations suggest that in 
some GCs the NS retention fraction could be as high as 
20% [T5]. Also, note that the tidal capture cross section 
used in [11] is more than an order of magnitude smaller 
than the GW capture cross section (discussed in the fol- 
lowing section) for compact objects, and using the latter 
would increase the rates by the same factor. In geometric 
units G = c = 1 (which, unless otherwise stated, we use 
throughout), tidal capture is estimated to occur in |llj 
for periapse values r p /M< 32, 25, and 13 for NS-NS, 
BH-NS, and BH-BH binaries, respectively. 

In [5U] NS-NS binary formation in GCs via exchange 
interactions was studied, giving a merger rate of ~ 2 yr _1 
Gpc~ 3 . A similar mechanism was explored in [21] for BH- 
NS systems; the results depend sensitively on the initial 
mass fraction of BHs, with more massive BHs leading 
to higher event rates. For example, models where the 
GC contained M = 35© BHs lead to advanced LIGO 
detection rates of 0.04-0.7 yr . Though in contrast to 
tidal or GW capture discussed in the previous paragraph, 
the mechanisms looked at in both these studies typically 
produce binaries with periods of 0.1 days or longer, and 
they will effectively circularize before entering the LIGO 
band. 

There is also the possibility that eccentric mergers 
could result from hierarchical triples through the Kozai 
mechanism. This has been suggested to occur in BH- 
BH mergers in GCs [HJ [221 [53] and CO mergers around 
supermassive BHs in galactic nuclei [53], as well as in 
coeval or dynamically formed BH-NS or NS-NS bina- 
ries [55]. Though the dynamics of these systems will be 
different from those studied here, they could be similar 
at late times. Efforts to understand this mechanism in 
the general-relativistic regime are ongoing (see e.g. [26]), 
and the event rates of these systems are not well known 
(though see [27]). 



B. Cross sections and high-eccentricity fractions 

We focus on mergers with initial periapse r p < 10M, 
where M is the total mass, making this study comple- 
mentary to previous studies [SJ [7]. As we show later, 
in this regime essentially all mergers occur with non- 
negligible eccentricity (e > 0.2). This is also the regime 
where strong-field effects such as black hole spin and 
zoom- whirl behavior can influence the dynamics. To es- 
timate the fraction of dynamical capture binaries that 
retain high eccentricity, we can use Newtonian dynam- 
ics with quadrupolar energy loss following [2"8"rl3"0"] . First, 
for a hyperbolic orbit with a small velocity at infinity 
»<1, the relationship between impact parameter b and 
r p is r p rj b 2 v 2 /2M. In other words, the cross section 
a oc b 2 scales linearly with r p . The maximum pericenter 
passage that leads to a bound system through gravita- 
tional radiation loss 1 is r Pim w (Sl-q) 2 / 7 v^ 4 / 1 M , where 
r\ = rnxm^jM 2 = q/(l + q) 2 is the symmetric mass ra- 
tio, with q the mass ratio. For a galactic nuclear cluster 
where v w 1000 km/s, between 20% and 30% of dynam- 
ical capture binaries (where the range is from q = 1 to 
q = 0.1) will have r p /M < 10; for a globular cluster with 
dwIO km/s, this drops to 1.5%-2.0%. 

Although we focus on those with small initial peri- 
apse, all dynamical capture binaries will have a repeated 
burst phase [7j- For a large fraction of expected binary 
masses the repeated bursts will be within the Advanced 
LIGO band. The burst frequency is r p 1 (r p /M)~ 1 / 2 ; 
the lowest frequency occurs at r p — r p m , which ranges 
from (1 — 100Hz)/Mio for q = 0.1 encounters in globu- 
lar clusters to q = 1 encounters in nuclear clusters, with 
Mm = M/IOM©. To estimate the percentage of systems 
that will end with a low-eccentricity inspiral phase, if the 
initial periapse is r p ^ , and we consider the repeated burst 
phase to end at a periapse of r p j with eccentricity e/, 

from [55] r p>i 0.57r pJ (l + e/ )e~ 12/19 [l + 0{e 2 )]. For 
example, if a binary with e/ < 0.1 by r p j = 10M can 
be considered to have a low eccentricity inspiral phase, 
then this corresponds to all systems with r Pt i > 27M. 
For nuclear clusters, this is between 20% and 40% for 
q = 1 — 0.1, while the corresponding range for globular 
clusters is 94% and 96%. 



C. Electromagnetic counterparts 

Binary NS or BH-NS mergers are thought to be pro- 
genitors for SGRBs, and may also source a number of 



As mentioned earlier, for the COs considered here, energy lost to 
tidal effects is much less than GW emission at these separations, 
so the latter process determines the cross section. Also, when a 
bound system is formed, the fraction that has a semimajor axis 
large enough to have the binary tidally unbound by a subsequent 
interaction with the surrounding cluster potential is insignificant. 
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other EM transients [3TJ [35] . Possibilities include opti- 
cal/UV emission on time scales of a day from radioac- 
tive decay of ejected material. (This depends on heavy 
element opacities. Recent work using more detailed cal- 
culations suggests the time scale may be up to a week 
with emission peaking in the IR [33 ) . Interaction of the 
outflow with surrounding matter can also produce radio 
emission on time scales of weeks to years [34J. And, for 
binary NS mergers, the strong shocks produced can emit 
in radio to x rays over a second to day time scales [35] . 

Simulations of eccentric BH-NS and NS-NS mergers 
have shown a rich variation in outcome with impact pa- 
rameter, with the possibility of large accretion disks as 
well as ejecta that could undergo the r process QTJ I36T - 
135] , There is also significant variability in observed 
SGRBs. It is not implausible that this may in part be 
due to a subclass of SGRBs associated with dynamical 
capture binaries. Though not conclusive, there is also 
observational evidence for multiple SGRB progenitors. 
Of SGRBs with identified host galaxies, ~ 25% have off- 
sets of > 15 kpc from their hosts [10], which would be 
consistent with kicked, primordially formed binary COs 
or with dynamically formed binaries in globular clus- 
ters. The latter may be preferred for the largest off- 
sets [U] , especially if primordial binary COs experience 
weak kicks [32]. X-ray afterglows suggest that differ- 
ent progenitors may be responsible for SGRBs with and 
without extended emission [43J ; simulations of dynamical 
capture binaries show that it is more common to get long 
tidal tails, which could lead to extended emission as the 
material falls back to the accretion disk. There is also a 
claim that a high-energy gamma-ray source observed in 
Terzan 5 may be the remnant of a binary CO merger- 
powered SGRB [33]; if true, this provides evidence that 
dense cluster environments can be significant sources of 
binary CO mergers. 

The time scales between close encounters in eccentric 
mergers may also explain observed delays between pre- 
cursors and SGRBs 45 . For example, NS crust cracking 
on a nonmerging close encounter could potentially cause 
flares that precede the merger by an interval ranging from 
milliseconds to possibly a few seconds [35]. 



D. Gravitational wave detectability 

Multimessenger exploration is of course an exciting 
possibility. Even a null GW detection provides astronom- 
ical information as it rules out compact object mergers 
as the source of an observed GRB, but only if the de- 
tectability of these types of signals is understood. Given 
the disparate nature of the waves from dynamical cap- 
ture vs coeval mergers, data analysis methods designed 
specifically for each are required for this kind of astron- 
omy. Methods to search for quasicircular inspiral (of rel- 
evance to the majority of coeval binaries, and a subset of 
dynamical capture binaries that form with a sufficiently 
large periapse to circularize before merger) have been 



the predominant focus of the GW community over the 
past decades [46]. Comparatively, there is a dearth of 
studies on the detectability of highly eccentric mergers 2 . 
In [45] the single burst from a parabolic close encounter 
was studied, while [6 included the additional signal pro- 
vided by subsequent bursts. This repeated burst phase 
was studied in [7] using 2.5 and 3.5 order post-Newtonian 
(PN) equations of motion. It was found that GWs from 
this phase may be detectable by Advanced LIGO out to 
200-300 Mpc for BH-NS binaries and 300-600 Mpc for 
BH-BH binaries. Since the PN approximations begin to 
break down close to merger, the evolution was only fol- 
lowed to r p = I OA/ . To model the last stages of merger 
requires numerical relativity (NR), and there have been 
a number of numerical studies of eccentric mergers |36l - 
1351 [45-54 . However, because of the computational ex- 
pense of these simulations, it is not possible with current 
computer resources to follow high-eccentricity binaries 
through multiple close encounters. The challenge is com- 
pounded by a large parameter space, including the im- 
pact parameter, mass ratio, BH spin, and NS equation of 
state. It is thus not reasonable to expect that brute-force 
numerical simulations will be able to provide templates 
before the Advanced LIGO era, even accounting for ex- 
pected increases in computer power. 



E. Outline of remainder of paper 

To begin to bridge the gap between large periapse PN 
solutions and late-time numerical solutions, we introduce 
a model for the inspiral, merger, and ringdown of dynam- 
ical capture compact binaries. This model is based on 
geodesic equations of motion in an effective Kerr space- 
time, combined with quadrupole radiation (Sec. II A) 
and a version of the Implicit Rotating Source (IRS) 
model [551 156| for the merger and ringdown parts of the 
GW signal (Sec. |II G\ . (Except for the IRS extension, 
and the comparable masses, our hybrid is reminiscent 
of the "kludge" introduced to study extreme mass ra- 
tio inspirals |57H60j , based in part on the "semirelativis- 
tic" approach of [ST].) We validate this model through 
a comparison to full numerical simulations in the strong- 
field regime (Sec. II B) and to the PN approximation for 
r p > 10 M (Sec.|HDf 

The waveforms we produce here are likely not accurate 
enough for optimal template-based detection of multiple- 
burst events. Indeed, creating improved accuracy wave- 
forms will probably require a different approach, for ex- 
ample, using the effective one body (EOB) formalism 
which has recently been extended to generic orbits [55] 
and calibrating it using full numerical simulations. How- 
ever, our waveforms capture the relevant features with 



2 Though see [47] for a recent study of the efficacy of quasicircular 
templates to detect lower eccentricity NS binaries. 
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sufficient faithfulness that we can use them to assess the 
efficacy of existing LIGO search strategies. We can also 
use our waveforms to investigate new search strategies 
that may be better suited to highly eccentric mergers. 
In Sec. IIIII we use this model to evaluate how well these 
GW signals could be seen with each generation of the 
LIGO detectors, varying impact parameter (equivalently 
r p ), total mass, and mass ratio. We use various analysis 
methods: matched filtering with the model templates, fil- 
tering with ringdown templates, and a burst search with 
sine-Gaussian templates. We also estimate how well a 
hypothetical search using incoherent stacking of bursts 
following [63 would perform. Though not as optimal as 
matched filtering, stacking is likely more robust to timing 
uncertainties in the burst sequence. We find that if cap- 
ture binaries do exist, in many cases their GW signals will 
be missed by single-burst or ringdown searches (and, as 
we argue, quasicircular templates), whereas these sources 
would be detectable with a full template or a stacked 
burst search. In particular, GRB051103 [53] had a mea- 
sured distance of 3.6 Mpc, and no coincident GW signal 
was found using traditional searches |651 166 ) . However, 
there is a sizable region of the parameter space of dynam- 
ical capture binaries that existing searches would have 
missed. The possibility that the GRB was preceded by 
an eccentric merger remains a viable possibility. 

In Sec. |IV| we make concluding remarks and comment 
on the direction of future work. 



II. WAVEFORM MODEL 

In this section we describe our model for high- 
eccentricity merger waveforms. We first look at the inspi- 
ral phase in Sec. |II A[ which can be considered a sequence 
of GW bursts, each generated at a periapse passage. In 
Scc. |II B| we compare the model expressions we use for the 
bursts to full numerical simulations. In Sec. Ill CI we dis- 
cuss the IRS model for the merger and ringdown phase, 



and in Sec. II D we present examples of the full signal, 
and make further comparisons to PN results for the in- 
spiral phase. 



A. Repeated burst phase 

Our objective here is to model the GW signal from 
an eccentric binary that passes through a series of close 
encounters prior to merger. To this end, we use a pre- 
scription based on the equations of motion of a geodesic 
in a Kerr spacetime, coupled with the quadrupole for- 
mula for gravitational radiation. We identify the mass 
and total angular momentum of the binary with the mass 
and spin parameters of the effective Kerr spacetime and 
the orbital angular momentum and energy with that of 
the geodesies. This approach has the advantage of re- 
producing the correct orbital dynamics in the Newtonian 
limit and general-relativistic test particle limit, while still 



incorporating strong-field phenomena such as pericenter 
precession, frame dragging, and the existence of unstable 
orbits and related zoom-whirl dynamics. For simplicity, 
in this first study we restrict our attention to equatorial 
orbits and, for the most part, nonspinning BHs (we com- 
pare the IRS model to a merger involving a spinning BH 
inSecliTC). 



The equations for an equatorial geodesic in a Kerr 
spacetime with mass M and dimensionless spin a can 
be written in first order form using Boyer-Lindquist co- 
ordinates as 



= Q, 



ER 2 - 2M 2 aL/r 



LQ + 2M 2 a/r 



Pr = 



R 2 

AQP r /r 2 , 

\n 2 (r 3 - M 3 a 2 ) + M(2Man - 1) 

r 2 Q 

+ ^ [M 2 a 2 -Mr] 



(1) 



where R a = r 2 + 2M 3 a 2 /r + M 2 a 2 , A = r 2 - 2Mr + 
M 2 a 2 , P r — r 2 r I (AQ), r is proper time, and the overdot 
indicates a derivative with respect to the coordinate time. 
Here E and L are the energy and angular momentum of 
the geodesic. 

In order to apply these equations to a binary system 
we go to the center-of-mass frame and let r be the separa- 
tion between the objects. Then we identify the geodesic 
parameters E and L with the reduced energy and angular 
momentum of the system and promote these quantities 
to time-dependent variables. To determine the amount 
of energy and angular momentum radiated away to grav- 
itational waves, we use the quadrupole formula 



E 
L 



[A... ... 



2/i 



<''-/, it 



jk 



(2) 



where \i is the reduced mass, is the reduced 

quadrupole moment, andZy and Xm are written in terms 
of the variables {r,<j),P r ,E,L} using We set M 

in ([!]) to the total mass (neglecting orbital energy con- 
tributions), and we set a = iiL/M 2 + clbh, where cibh 
is the net spin of any BHs (though again for this study 
we focus on nonspinning BHs, where asn = 0). The 
use of an effective spinning BH spacetime based on total 
angular momentum is motivated by [49j . where it was 
found that the properties of zoom-whirl-like dynamics 
exhibited in equal mass mergers in full numerical rela- 
tivity are better approximated by geodesies on the effec- 
tive Kerr spacetime than Schwarzschild spacetime, and it 
differs from the EOB approach which uses deformations 
of the Schwarzschild metric for the merger of nonspin- 
ning objects [57] • We note that when the orbital angular 
momentum is large we will have a > 1. However, this 



5 



will occur only when the separation r is also large, so 
general-relativistic effects are small, and no unusual be- 
havior arises from exceeding the Kerr limit. We numeri- 
cally integrate the coupled set of equations (I]) and (j2J). 

The remaining element is to calculate the observed 
gravitational radiation, which will depend on the intrin- 
sic source parameters (i.e. the mass, mass ratio, eccen- 
tricity, and initial periapse distance), and will also vary 
with sky location and relative orientation of the source to 
the detector. At linear order and in the transverse trace- 
less gauge, the complex gravitational wave strain h opt a 
distance d from an optimally oriented source is simply 
related to changes in the quadrupole moment through 

^Opt = + th Opt = l(l xx + l Xy „) . (3) 

For general orientations, the emitted strain can be rep- 
resented through a mode decomposition as 
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TABLE I: Fit parameters for close-encounter GWs. 

a Approximate initial parameters of the geodesic model based 

on the initial orbital energy and angular momentum of the 

simulation. 

b Initial parameters of the geodesic model that best fit the 
simulation data. 

c Amplitude enhancement applied to the waveform from the 
best-fit geodesic model. 

d Overlap between simulation and the best-fit geodesic model. 



E 



E 

1=1 m=-t 



h em (t,d)- 2 Y £m (9,(b), (4) 



where - 2 Yi m are the spherical harmonics of spin weight 
—2 [68] . and and 4> are the polar and azimuthal angles of 
orientation, respectively. For the comparable mass, non- 
spinning systems that we are primarily interested in, the 
quadrupole (i.e., I = 2, m = +2) component dominates 
the strain, so that 

h « h 22 (t, d) - 2 Y 22 {0, 4>) + h 2 -2(t, d) _ 2 F 2 _ 2 (0, <f>). (5) 

This completes the approach for calculating the source 
waveform that reaches a detector. In a later section we 
will include the sensitivity of the detector in the analysis. 



Comparison to fully general-relativistic 
numerical simulations 



To provide some validation for this model we compare 
several waveforms of single high-eccentricity fly-by en- 
counters from full general-relativistic numerical simula- 
tions to those obtained from the geodesic equation with 
the quadrupole formula. The simulations include a 4:1 
mass ratio BH-NS system [37] , an equal mass NS-NS sys- 
tem [55], and an equal mass BH-BH system. The NR 
simulations were all performed using the code described 
in |52| . 

In Fig. [T] we show several such examples from NR simu- 
lations of the 4:1 BH-NS system alongside corresponding 
waveforms from our model with best-fit parameters. The 
peak amplitude of the geodesic is scaled to be the same 
as in the simulations. The fit is performed by finding the 
initial orbital parameters that maximize the phase over- 
lap between the waveforms (see e.g. |70j). In this regime 
the match between the waveforms is most sensitive to r p 
as opposed to e. As can be seen, the fly-by waveforms 
from our model provide a good match to those from sim- 
ulations. Even close to the effective innermost stable 



orbit (ISO) for the BH-NS system (the bottom panels 
of Fig. [I]), where the system begins to show evidence of 
whirling behavior, our model is able to approximately 
capture the shape of the waveform. 

In Table |T] we give the fit parameters, amplitude en- 
hancement, and overlap. We also show the approximate 
initial orbital parameters (r p and e) of the simulation ob- 
tained by equating a Newtonian estimate of the reduced 
orbital energy and angular momentum at the beginning 
of the simulation with the E and L parameters of the 
geodesic model described above. (Note, this is different 
from the Newtonian values for r p and e used in |37U38| .) 
For most of the BH-NS systems in Table [I] we can see 
that the enhancement required to match the amplitude 
of our model to the simulation results is ~ 4%-ll%. This 
is presumably due to aspects not captured by this simple 
model, such as finite-size effects, as well as truncation 
error from the simulations. The one case where the am- 
plitude of the simulation waveform was below the model 
result was a simulation with strong whirling behavior 
(bottom-right panel of Fig. [I]) where the NS had large 
/-mode oscillations excited as described in [37] . 

We also compare the geodesic model with an equal 
mass BH-BH and an equal mass NS-NS system as shown 
in Fig. [2] Although one would expect a geodesic approx- 
imation to be most accurate in the limit that one mass is 
much larger than the other, it still provides good fits for 
equal masses. This model, however, does not attempt to 
include finite-size effects (such as the /-mode excitation 
visible in the later part of the bottom of Fig. [2]), which 
would be required to address questions related to mea- 
suring the NS equation of state from such GW signals. 
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FIG. 1: Comparison of the £ = 2, m — 2 component of h for fly-by waveforms from 4:1 BH-NS simulations (solid line) and 
our model with best-fit periapse distance (dashed line). The approximate effective geodesic orbital parameters of the simulated 
system (left to right, top to bottom) are (r p , e) = (8.3, 1.0), (8.0, 0.8), (5.6, 1.0), and (5.0, 1.0). The fit parameters are given in 
Table H] 



C. Merger model 

After a binary has evolved through some number of 
close encounters, it will merge. In order to include the 
waveforms resulting from merger, we supplement the 
model outlined in Sec. Ill Al with a version of the IRS 
model [55J for the merger and ringdown part of the 
GW signal. Note that the IRS assumes the waveform is 
circularly polarized. This is not strictly valid for the 
complete merger-ringdown phase of eccentric binaries, 
though as we show below, it does provide a reasonably 
good approximation to results from numerical simula- 
tions. As with other aspects of our waveform model, this 
assumption could be refined in the future, but it is ad- 
equate for the purpose of testing the efficacy of existing 
search strategies for detecting eccentric binaries. 

In particular, we model the phase evolution to asymp- 
totically approach the least damped quasinormal mode 
frequency of the final BH, wqnm> via 



w(t) =wqnm(1 -/) (6) 



where 

/=^(l + -) 1+K fl-(l + -e- 2t/b )-' i V (7) 

2 K \ K J 

Here b = 2Q/wqnm is determined by the quality factor 
and frequency of the final BH, and k and c are free pa- 
rameters of the model. The amplitude is given, up to an 
overall factor Aq, by 



w(*) ^l + a(/ 2 -/ 4 ) J 

where / = df/dt, and a is a free parameter. We find 
that a = 72.3/Q 2 provides a reasonably good fit to our 
numerical simulations. 

In Fig. [3] we show a comparison between simulation 
results of BH-NS mergers and the best-match IRS model 
waveforms, where we let k and c be fitting parameters. 
In Fig. [4] we show the same thing for equal mass NS-NS 
and BH-BH mergers. This simple model will not capture 
disruption or other matter effects, and best-fit values for 
k and c will have some dependence on the parameters 
of the binary, such as the impact parameter preceding 
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FIG. 2: Comparison of the £ — 2, m = 2 component of 
h for fly-by waveforms from equal mass BH-BH (top panel) 
and NS-NS (bottom pannel) simulations (solid lines) and our 
model (dotted lines) with best-fit parameters. The approx- 
imate effective geodesic orbital parameters of the simulated 
systems are (r p , e) = (8.7, 1.0) for both cases. The fit param- 
eters are given in Table [I] The feature in the waveform after 
the peak in the NS-NS simulation is from /-mode excitation 
that occurs during the close encounter. 



merger. However, when studying signal detectability we 
fix k — 0.64 and c = 0.26, which empirically provides 
reasonably good fits to a large number of simulated wave- 
forms, and therefore provides an adequate representation 
of a generic eccentric merger. We attach the IRS part of 
the waveform to the model from Sec. |II A] when the sep- 
aration reaches the light ring of the effective Kerr space- 
time. 



D. Model properties and comparison to 
post-Newtonian 

Combining the inspiral and merger models allows us 
to generate complete waveforms for dynamical capture 
binaries. In Fig. [5] we show one such example for a 4:1 
mass ratio system with initial orbital parameters corre- 
sponding to r p = 8 M and e = 1. The waveform shows 
the decreasing time interval between bursts from close 
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FIG. 3: Comparison of the merger GW strain from a 4:1 BH- 
NS simulation (solid lines) and the IRS model (dotted lines) 
with best-fit parameters. The top panel shows a case where 
the initial BH was nonspinning. The bottom panel shows a 
case with obh = 0.5, which results in more whirling behavior 
and tidal disruption of the NS. The best-fit parameters in Q 
are (k,c) = (0.66,0.28) (top) and (0.46,0.18) (bottom), and 
the matches are 0.98 and 0.96, respectively. The match is 
weighted based on the "whitened" waveforms as described in 
Sec. |III| assuming a total mass of 10 Mq. 

encounters as r p and e decrease due to gravitational radi- 
ation. The number and timing of the bursts is a sensitive 
function of the amount of energy and angular momentum 
radiated in each close encounter. In Fig. [6] we show how 
r p and e evolve according to this model for some example 
binaries. It can be seen that the binaries considered here, 
which begin on parabolic orbits with r p <10M, still have 
non-negligible eccentricity all the way to merger. 

We can also compare this model to that given by the 
2.5 and 3.5 order PN approximation as used in [TJ. In 
Fig. [7] we show how the difference in the energy and an- 
gular momentum radiated away in a close encounter for 
2.5 or 3.5 PN relative to our model changes with the ini- 
tial impact parameter. The geodesic model predicts less 
energy and momentum loss than 2.5 PN but more than 
the 3.5 PN. At large impact parameters the three dif- 
ferent models converge. At smaller impact parameters 
the 2.5 and 3.5 PN approximations begin to diverge. As 
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FIG. 4: Comparison of merger waveforms from an equal 
mass BH-BH simulation (top panel) and NS-NS simulation 
that forms a BH (bottom panel) with the IRS model (dotted 
lines) with best-fit parameters. The best-fit parameters in |7|) 
are (k,c) = (0.31,0.03) (top) and (0.36,0.19) (bottom), and 
the matches are 0.98 and 0.97, respectively. The match is 
weighted based on the "whitened" waveforms as described in 
Sec. |III| assuming a total mass of 20 Mq and 2.8 Mq for the 
BH-BH and NS-NS binaries, respectively. 
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FIG. 5: GW strain generated with our model and initial 
conditions r p = 8 M and e = 1. The top panel shows the 
entire waveform, while the bottom panel shows a zoomed-in 
view of the end of the waveform. 



III. DETECTABILITY 
A. Detector modeling 



shown in [7T] , the PN approximation fails to converge (or 
even to provide physically sensible results in the case of 
3.5 PN) for r p < 10 M. 

The gravitational wave model we have outlined in this 
section is relatively simple and could be improved upon 
by, for example, adding more sophisticated conservative 
dynamics, including finite-size effects for NSs, as well as 
going beyond the quadrupole approximation in determin- 
ing gravitational radiation. However, given the decent 
match between this model and the full numerical simu- 
lations, as well as its consistency with the PN approxi- 
mation as described above, it can be used to investigate 
issues of detectability, as we do in the next section. 



Having developed a model for the gravitational wave- 
forms emitted by high-eccentricity binaries, we can now 
assess the detectability of these signals for different 
source parameters and detectors. The measured strain 
h is given by 

h = [Fh] = F+h+ + F x h x , (9) 

where F = F + — iF x is the sky-dependent detector re- 
sponse. The signal-to-noise ratio (SNR) p using a per- 
fectly matched filter is given by 

p 2 = (h\h), (10) 

where (-|-) denotes a noise- weighted inner product given 
by 

+ (11) 

o 
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FIG. 6: Evolution of orbital parameters for a 4:1 (top panel) 
and an equal mass ratio (bottom panel) binary. The effective 
eccentricity is calculated from successive apoapse and peri- 
apse distances as e = (r a — r p )/(r a + r p ). 
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FIG. 7: Relative difference in energy and angular momentum 
lost in a close encounter for the 2.5 and 3.5 PN approxima- 
tions versus our model. For this comparison the orbits are 
chosen to initially have zero energy and the same value of an- 
gular momentum which, for a geodesic with the same initial 
conditions, corresponds to the value of r p indicated on the x 
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FIG. 8: The characteristic strain h c is shown for the ini- 
tial (thin dash-dotted line), Enhanced (dash-dotted line), and 
Advanced LIGO (solid line) detectors, as well as for two ex- 
ample signals at Dl = 1 Gpc. The first signal corresponds 
to an orientation-averaged source with M = 100 Mq, 9=1, 
and r p = 5M (dashed line), and the second signal is from a 
source with M = 10 M Q , q = 0.1, and r p = 10M (dotted line). 
Both signal spectra are smoothed to diminish fluctuations and 
make the trend more clear. The system with q — 0.1 has lit- 
tle contribution from the merger, so the repeated burst phase 
dominates the spectra, with h c %f ', whereas the q — 1 system 
signal comes largely from the merger, where h c m constant 
over a small band of frequencies. 



where S n (/) is the power spectral density of the detector 
noise, and h denotes the Fourier transform of the orig- 
inal h time series. Because we limit our model to the 
quadrupole component of the signal, and we focus on 
detectors (like LIGO) for which the gravitational wave- 
length is much longer than the detector's armlength, we 
can trivially relate the SNR of an optimally oriented 
and located source to the SNR of an orientation- and 
sky-location-averaged source. For such detectors, the re- 
sponse function to the two waveform polarizations, F + 
and F x , is simply the root-mean-squared (rms) average 

over the sky location and polarization angles, yj (F 2 x ) = 

y/l/5 [75]. Likewise, the rms average over source orien- 

/1/5, so that d (p 



tations is 



/5- 



'2 IT 2' v ~/ — V \P I Popt/ 

We can further define the characteristic strain h c for both 



the signal and detector noise. Given Eq. ( 11 ) and the typ- 



ical practice of plotting sensitivity curves logarithmically, 
it is useful to define h r = 



? )fh op t for signals 



and h c = yj fS n /{F 2 ) for detector noise, so that both sig- 
nal and noise are characterized as a dimensionless strain, 
and the ratio of signal-/i c to noise-/i c is the square root 
of the integrand for p 2 when integrated over logarithmic 
frequency intervals df / f. We show this characterization 
of signal and noise in Fig. [8] 

For assessing the relative contribution of different 
waveform segments to the SNR, it is often convenient 
to work in the time domain by constructing "whitened" 
waveforms 73] , which weight the amplitude of the wave- 
form as a function of frequency to account for the pres- 
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FIG. 9: Whitened waveforms for a 10 Mq (top panel) and 
a 100 M0 (bottom panel) binary with initial e = 1 and r v = 
5M, along with the (whitened) best-fit template among sine- 
Gaussian and ringdown templates. 



ence of noise in the detector, 

+00 

h'= I df- 



-i2wft 



(12) 



With these whitened vectors, the noise-weighted inner 
product (11 1 can be reexpressed in the time domain: 



{h!\h 2 ) = / dth'*(t)h' 2 (t). 



(13) 



Figure [9] shows portions of the whitened waveform for 
two example cases with the same mass ratio and initial 
r p and e, but different masses. The upper panel shows 
the burst with the largest SNR contribution for a source 
with total mass M — 10 M , while the lower panel shows 
the loudest burst for M = 100 M . The different masses 
change the frequency of the signal, so different bursts are 
emphasized by the detector sensitivity; in particular, for 
larger masses the final burst and merger are emphasized. 
We also show best fits for two types of templates that are 
described below. 



B. Templates and detection strategies 

While quasicircular sources are searched for using 
matched filtering, eccentric systems are far more sus- 
ceptible to modeling errors in the relative timing and 
phase of signal bursts, which is why we focus our atten- 
tion on alternative approaches to detection. For example, 
a small modeling error in the energy lost during a par- 
ticular periapse passage SE will induce a timing error 



in the arrival time of the subsequent burst 8T given by 
ST oc SE(1 — e)~ 5 / 2 . Therefore, dynamical capture bi- 
naries are far more challenging to model with sufficient 
accuracy to apply matched filtering due to their large 
eccentricities. 

We assess detection prospects of GWs from capture 
binaries for two currently used templates, sine-Gaussian 
(SG) and ringdown (RD) templates, as well as an ideal- 
ization of a third strategy based on combining an excess 
power search with stacking. The SG and RD both take 
the form 



h = A exp 



t-t c 



+ iuj(t — t Q ) — i(f 



(14) 



where 7 = 1 and t>t D for the RD templates, and 7 = 2 
with -co < t < 00 for the SG templates. Here A is the 
overall amplitude, t a and <p are the time and phase of 
the template's amplitude peak, r sets the e-folding of the 
amplitude, and uj is the constant frequency. 

In addition to assessing the performance of two burst 
templates, we calculate a rough approximation of the po- 
tential performance for an excess power search that ac- 
cumulates power from the entire signal [63] . which we 
will call a power-stacking search. Here, the data would 
be transformed to a time-frequency (TF) tiling using a 
basis suitable for capturing individual bursts within a 
tile, and then power from different tiles corresponding to 
bursts, as informed by our model, would be combined. 
Whereas in most existing TF searches an individual ele- 
ment must have enough power to exceed some threshold, 
with that threshold being large enough to avoid many 
false alarms, the approach we describe does not require 
that the signal be detectable in any single TF element. 
In the case of a monochromatic signal, the SNR from 
optimal filtering will accumulate with the number of cy- 
cles N as y/N, while the excess power in stacking TF 
elements (constructed using any basis) overlapping the 
given frequency will accumulate as N 1 / 4 . The signals 
from high-eccentricity binaries are not monochromatic, 
but given the typically large number of bursts occurring 
in band, and the relative flatness of both the source spec- 
trum and the detector sensitivity across its most sensitive 
band, we expect the aforementioned scalings to hold ap- 
proximately for realistic signals. 

This search would be very similar to the stacked search 
proposed for combining potential GW counterparts to 
observed electromagnetic signals from soft gamma re- 
peaters |63j . There, TF elements were aligned in time 
based on the observed bursts, and they demonstrated 
the iV 1 / 4 SNR scaling when adding power for identical 
injected signals. Since we do not have a separate obser- 
vational trigger, our proposed search would sum power 
along elements overlapping bursts as indicated by our 
waveform model. We leave it to future work to fully 
investigate this, though here we assume we can achieve 
the N 1 / 4 SNR scaling, and thus can estimate the per- 
formance of a power-stacking search by noting that opti- 
mal filtering should outperform power stacking by N 1 / 4 . 
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Hence, we can approximate the effective excess power 
SNR as pep ~ N~ 1 / 4 p. This simpie estimate will consti- 
tute our third search technique in our subsequent analy- 
sis. 

We do not employ quasicircular (QC) templates, al- 
though they have thus far been the only tool employed 
to search for long-lived signals. QC templates will gener- 
ically fail to match the performance of any of the above 
methods for the repeated burst phase of eccentric sources 
for the following reasons. First, during the long intervals 
between eccentric bursts a QC template will still be in- 
tegrating power from the data, which is predominantly 
noise. Specifically, the ratio of the characteristic time 
scale of an eccentric burst to the period between bursts 
is roughly 



T~GW 

T 



(1-e) 



3/2 



(15) 



In other words, there will be ~ (1 — e)~ 3 / 2 additional 
cycles between bursts in a QC signal with the same peri- 
apse. Moreover, even if the QC template is phase aligned 
to a particular burst, since the time between bursts is 
much larger than the GW period, the rest of the tem- 
plate will effectively have random phase alignment with 
other bursts in the sequence and, on average, no addi- 
tional SNR will be acquired. To summarize, typically 
the best-matched QC template will only integrate signal 
about the loudest burst, but even so, the performance 
will not be as good as a single-burst search due to the 
larger integrated noise accumulated over the period of the 
QC template (expect for the higher mass systems where 
only the final merger/ringdown signal is in band). 



C. Results 

We calculate two useful quantities related to the SNR: 
the detectability horizon and the detection probability. 
Since h oc DJ 1 , where Dl is the luminosity distance, 



we can use ( 10 1 to calculate the distance (which we call 



the detection horizon) at which a sky- and orientation- 
averaged source could be observed with a SNR of 8 using 
optimal filtering. The detection probability for a given 
strategy is simply the ratio of the volume in which the 
strategy could detect a source with some SNR to the 
volume in which the source could be seen with the same 
SNR using optimal filtering. In the remainder of this 
section, we will calculate these quantities for various cases 
of interest. We consider the following configurations: 

• three detector sensitivities, corresponding to initial, 
Enhanced, and Advanced LIGO; 

• three detection strategies, including SGs, RDs, and 
power stacking, and how they compare to optimal 
filtering: 

• three intrinsic system parameters: 



— total system mass M, ranging from 1 M Q - 
2000 M ; 

— mass ratio q of the binary components, rang- 
ing from 0.01-1; 

— initial r p , ranging from 5M-10M (with initial 
e = 1; we exclude r p < 5M simply because in 
most cases it is a direct collision qualitatively 
similar to r p — 5, and see [7] for a study of 
r p > 10M). 

In Fig. [lOj we show contours of constant horizon dis- 
tance as a function of q and r p for initial LIGO, assum- 
ing optimal filtering, SG templates, and power stacking. 
Two contours of note, at 0.77 and 3.6 Mpc, correspond 
to the distances of GRB070201 [7J and GRB051103 [55, 
respectively. These were two nearby gamma-ray bursts 
observed by Swift during the S5 initial LIGO run, while 
two interferometers were actively collecting data at or 
near initial LIGO's design sensitivity. However, no sig- 
nal was found in the LIGO data using the methods ap- 
plied (specifically, various burst and quasicircular inspiral 
templates) for these GRBs, nor for any of the 137 GRBs 
(35 with measured redshifts) that occurred while initial 
LIGO was taking science data during its S5 run at or 
near design sensitivity 
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we re- 



[651 E6J. Thus in Fig. _ 
strict the mass ratio to the range 0.1-1, with one of the 
masses fixed at 1.35 M , to focus on systems including 
a neutron star that are expected to generate GRBs. In 
the case of a dynamical capture binary source at 0.77 
Mpc, the signal is sufficiently loud that even suboptimal 
searches like SG templates would detect them. However, 
for a source at 3.6 Mpc, whereas an optimal filter would 
have detected a signal from a large region of the parame- 
ter space, including all cases with q < 0.5 or r p > 7.5 M, 
and power stacking would recover signals with q < 0.4, 
SG templates are far less effective, and would only re- 
cover a small sliver of parameter space with q < 0.2. 
This suggests the possibility that the searches applied to 
the LIGO data would not have found the gravitational 
wave counterpart to GRB051103 if it was in the form of 
a dynamical capture binary. Furthermore, across the full 
parameter space explored, the difference in performance 
among these three searches is substantial, with optimal 
filtering detecting sources as far as Dl = 50-100 Mpc, 
while power stacking only reaches Dl ss 30 Mpc, and SG 
templates only reach Dl ss 15 Mpc. 

Figures [TT] and [12] show, for Enhanced and Advanced 
LIGO respectively, contours of detection horizon as a 
function of mass and mass ratio at a fixed r p — 6 M 
using an optimal filter and SG and RD templates. The 
primary difference in both cases is the degradation of per- 
formance for higher mass ratios (smaller q), with the SG 
performing as well as or better than the RD templates 
across much of the parameter space, with the exception 
of comparable mass ratios, where the ringdown signal 
is most emphasized. For each search, Enhanced LIGO 
could detect an equal mass binary with M = 100 M Q 
out to Dl — 1 Gpc, and Advanced LIGO will see the 



12 



0.5 0.77 1 2 3.6 5 10 20 50 100 



1 3 10 30 100 300 1000 



D L [Mpc] 



0.6 



0.2 




5 6 


7 8 

r P [M] 




9 


10 


0.5 0.77 1 2 


3.6 5 10 


20 


50 


100 





D L [Mpc] 




a-10 



cr10 



'P [M] 




1000 



1000 



FIG. 10: Contours of horizon distance (p — 8) as a function 
of mass ratio q and pericenter separation r p for initial LIGO 
using an optimal filter (top panel), sine-Gaussian templates 
(middle panel), and an estimate of a power-stacking search 
(bottom panel) as described in the text. We fix one com- 
ponent to be a 1.35 Mq neutron star and change the total 
mass with mass ratio accordingly. We include a contour at 
Dl = 0.77 Mpc and another at 3.6 Mpc, to show the region 
of parameter space where existing LIGO searches would not 
have seen a gravitational wave counterpart to GRB070201 [74] 
and GRB051103 [64], respectively. 



same sources beyond 10 Gpc. 

The relative performance of SG and RD is further 



demonstrated in Figs. 13 and 14 which show the de- 
tection probabilities of each template (equivalently, the 
ratio of the detectable volume using the templates to the 
volume using optimal filtering). SG templates perform 



FIG. 11: Contours of horizon distance as a function of 
rest mass M and mass ratio q for Enhanced LIGO using an 
optimal filter (top panel), sine-Gaussian templates (middle 
panel), and ringdown templates (bottom panel) for an initial 
pericenter separation of r v = 6M. 



best for M ~ 200 M systems using Enhanced LIGO and 
M 1000 M Q systems using Advanced LIGO, largely in- 
dependent of the mass ratio. Interestingly, RD templates 
perform best for comparable mass binaries regardless of 
total mass for Enhanced LIGO, whereas no such clear 
general behavior is observed for Advanced LIGO. This 
can be understood because Enhanced LIGO always has 
fewer cycles in band than Advanced LIGO, so that the 
merger-ringdown constitutes a larger fraction of the total 
SNR, with that fraction further enhanced for compara- 
ble masses (since p oc r\ — mim^/M for inspirals, but 
p oc ^frj for ringdowns [75]) . Advanced LIGO shows no 
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FIG. 12: Contours of horizon distance as a function of 
rest mass M and mass ratio q for Advanced LIGO using an 
optimal filter (top panel), sine-Gaussian templates (middle 
panel), and ringdown templates (bottom panel) for an initial 
pericenter separation of r v = 6M. 



such behavior because the number of inspiral cycles is 
so large that the merger-ringdown rarely dominates the 
total SNR. 

Figures [15] and [16] again show contours of horizon dis- 
tance for Enhanced and Advanced LIGO, but as a func- 
tion of total mass and initial pericenter distance at fixed 
q = 1. Both q and M rather strongly affect the de- 
tectability of sources over the range of masses considered. 
r p moderately affects the detectability for lower mass sys- 
tems (M < 20 M Q ), though very little for higher mass 
systems (which is expected since the number of bursts 
varies significantly with r p in the range 5M < r p < 10M, 
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FIG. 13: Contours of detection probability p = V/V m ax as a 
function of rest mass M and mass ratio q for Enhanced LIGO 
for a source inside the optimal filtering distance horizon, us- 
ing sine-Gaussian (top panel) and ringdown (bottom panel) 
templates for an initial pericenter separation of r v = %M . 



but as the mass increases fewer of the initial bursts are 
in band). SG templates outperform RD templates for all 
but the extremely high-mass systems, and a small region 
of extremely low-mass systems with very small r p , t hat 
merge after 0(1) orbit. This is also clear in Figs. 
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and 

|18[ which show the corresponding detection probabilities. 
In addition to SG and RD templates, Figs. [16] and [T8] 
show the relative performance of a stacked power search, 
which readily outperforms burst template searches for 
the full range of parameters. Since this is the case for 
q = 1, it will apply more so for cases with q < 1, as 
they experience more cycles, so that we can conclude that 
a power-stacking search will always outperform a burst 
search and is likely to be the optimal search approach in 
the absence of a matched-filter bank. SG and RD tem- 
plates perform best for M in the range 100 M Q -200 M 
for both Enhanced and Advanced LIGO, with the range 
of horizon distances being the same as in the M — q plots. 
This is as expected given that setting q — 1 maximizes 
the signal power at fixed M and r p . 

As shown in Fig. [18] all three search methods approach 
optimal-filter performance for large masses, M>500 M , 
since all three methods benefit from having the SNR con- 
centrated in a small number of cycles. However, for lower 
masses and therefore a larger number of in-band cycles, 
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FIG. 14: Contours of detection probability p = V/V ma x as a 
function of rest mass M and mass ratio q for Advanced LIGO 
for a source inside the optimal filtering distance horizon, us- 
ing sine-Gaussian (top panel) and ringdown (bottom panel) 
templates for an initial pericenter separation of r v = 6M. 



the SG and RD template performances degrade much 
more rapidly than power stacking. Whereas SG and RD 
templates reach detection probabilities as low as 0.01%, 
power stacking remains above 10% for the full parameter 
space considered. Since our power-stacking estimate is 
an idealization, uncertainties in the timing and frequen- 
cies of eccentric bursts may degrade the performance of 
a true TF power-stacked search. On the other hand, the 
results of [S3] suggest this method is rather robust to tim- 
ing uncertainties that are smaller than the characteristic 
time of each burst . 
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FIG. 15: Contours of horizon distance as a function of rest 
mass M and pericenter separation r p for Enhanced LIGO 
using an optimal filter (top panel), sine-Gaussian templates 
(middle panel), and ringdown templates (bottom panel). The 
mass ratio is q = 1. 



IV. CONCLUSIONS 

We have developed a novel waveform model for eccen- 
tric binary gravitational waveforms which can be applied 
for r p <10M, where conventional post-Newtonian wave- 
forms fail. Such binaries may form through dynamical 
capture in dense stellar environments. Our model is not 
sufficiently accurate to generate a matched-filter bank, 
and doing so will be very challenging for large eccentrici- 
ties. However, the model is adequate to supply mock sig- 
nals to explore the performance of existing LIGO searches 
in detecting highly eccentric biliary systems. Of existing 



search strategies, the ringdown and burst searches are 
best adapted to these systems. However, we find that 
a large fraction of the parameter space, where we in- 
cluded the impact parameter 5 < r p < 10M (see [7] 
for a complementary study of r p > 10M), total mass 
M e [1, 2000] M and mass ratio q £ [0.01, 1], has a sig- 
nificantly smaller horizon distance than what is, in princi- 
ple, achievable with a matched-filter search. This implies 
that a corresponding volume of sources could have been 
missed in prior searches and may be missed in future 
searches if better adapted strategies are not employed. 
Though it may be impractical to construct templates 
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in the near future (via numerical or analytical methods) 
for these systems that are accurate enough for optimal 
searches, a refinement of the waveform model presented 
here should be adequate for informing a power-stacking 
search. This method has the potential to increase SNR 
by « N 1 / 4 for an iV-burst event compared to a single- 
burst search. Though less than the effective iV 1 / 2 scaling 
of a full template search, this would still be a significant 
improvement. Note also that even for systems with larger 
impact parameters that do evolve to an essentially qua- 
sicircular inspiral following the burst phase, for most ex- 
pected binary parameters the burst phase will be within 
the band of ground-based detectors. Thus, the quasi- 
circular inspiral phase will be truncated compared to a 
primordial quasicircular inspiral, and though such a sys- 
tem may still be detectable with a quasicircular template, 
it would of course be misidentified, and a bias would be 
introduced in the estimation of the binary parameters. 

For future work, we intend to implement a power-stack 
search using this waveform model to fully explore the ef- 
ficacy of this method and its (in) sensitivity to timing 
errors, as well as continue to refine the model to in- 
clude (for example) spin precession and finite body ef- 
fects for neutron stars. We mentioned that the standard 
PN equations are ill suited to studying the late stages 
of mergers, in particular for high-eccentricity binaries, 
motivating our development of the effective Kerr with 
radiation-reaction model described here. However, the 
EOB approach [67] is an alternative to the PN expansion 
that is well behaved all the way to merger for quasicircu- 
lar orbits. This approach has recently been extended to 
generic orbits [05], and it will be interesting to explore 
EOB as the basis for a repeated burst waveform model. 
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FIG. 16: Contours of horizon distance as a function of rest 
mass M and pericenter separation r p for Advanced LIGO 
using, from top to bottom, an optimal filter, sine-Gaussian 
templates, ringdown templates, and a power-stacking search. 
The mass ratio is q = 1. 
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FIG. 17: Contours of detection probability as a function of 
rest mass M and pericenter separation r p for Enhanced LIGO 
for a source inside the optimal filtering distance horizon, us- 
ing sine-Gaussian (top panel) and ringdown (bottom panel) 
templates. The mass ratio is q — 1. 
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FIG. 18: Contours of detection probability as a function 
of rest mass M and initial pericenter separation r p for Ad- 
vanced LIGO for a source inside the optimal-filtering hori- 
zon distance, using sine-Gaussian (top panel) and ringdown 
(middle panel) templates and a power-stacking search (bot- 
tom panel). The mass ratio is q = 1. Note the different scale 
in the bottom figure. 
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